Domain wall dynamics in magnetic nanotubes driven by an external magnetic field
Li Zai-Dong1, 2, †, Hu Yue-Chuan1, He Peng-Bin3, Sun Lin-Lin4, ‡
Department of Applied Physics, Hebei University of Technology, Tianjin 300401, China
Key Laboratory of Electronic Materials and Devices of Tianjin, School of Electronics and Information Engineering, Hebei University of Technology, Tianjin 300401, China
School of Physics and Electronics, Hunan University, Changsha 410082, China
School of Public Management, Beihang University, Beijing 100191, China

 

† Corresponding author. E-mail: lizd@hebut.edu.cn vickysunlin@126.com

Project supported by the National Natural Science Foundation of China (Grant No. 61774001), the National Social Science Foundation of China (Grant No. 17BJY103), the Key Project of Scientific and Technological Research in Hebei Province, China (Grant No. ZD2015133), and the Construction Project of Graduate Demonstration Course in Hebei Province, China (Grant No. 94/220079). Peng-Bin He was supported by the Natural Science Foundation of Hunan Province, China (Grant No. 2017JJ2045).

Abstract

We use the Landau–Lifshitz–Gilbert equation to investigate field-driven domain wall propagation in magnetic nanotubes. We find that the distortion is maximum as the time becomes infinite and the exact rigid-body solutions are obtained analytically. We also find that the velocity increases with increasing the ratio of inner radius and outer radius. That is to say, we can accelerate domain wall motion not only by increasing the magnetic field, but also by reducing the thickness of the nanotubes.

1. Introduction

Owing to the tubular structure, the magnetic nanotubes have different properties compared with other nanomaterials. For magnetic nanotubes, magnetism is one of the main features and its magnetic properties can be used to record information, information transmission,[13] and realize logic devices.[4,5] In recent years researchers have synthesized many types of nanotubes[69] through continuous exploration, and often use the inside radius, outside radius, and length to describe the nanotubes. In the application of biology, both inside and outside nanotubes have strong surface activity, which can be used in other substances, adsorption in drug load, targeted therapy, and cell culture.[10] In the application of chemistry, the magnetic nanotube has a catalytic performance which is far greater than other nanoparticles. Hence, the magnetic nanotubes in catalytic chemistry have important research significance.[11]

In the research of nanomaterial magnetic problems, the domain wall (DW)[1215] in the nanostructures plays a very important role. Most of the previous studies on DW propagation in magnetic nanostructures focused on flat strips.[1619] A great deal of research has already been done, but most are related to one- and two-dimensional nanostructures, while the research on three dimensional nanotubes is lacking. In general, the motion of a DW can be driven by a magnetic field[1923] or spin-polarized current.[2427] On the flat strips, the well-known Walker breakdown of DWs[28] presents a major obstacle. Above a critical velocity, when the Walker breakdown occurs, the original DW structure collapses and it results in an abrupt drop of DW velocity. Still under discussion is whether there is an alternative to the Walker breakdown on the magnetic nanotubes. Analytic models by Landeros and coworkers predicted a decrease in velocity after a certain threshold.[29,30] Simulations by Yan et al. claimed the suppression of the breakdown, at least in certain geometric dimensions.[31] If this is the case, the domain wall velocities could reach the phase velocity regime, which in turn would lead to Cherenkov-like spin wave emission.[32] Magnetic nanotubes belong to a new class of interesting materials where the surface curvatures strongly influence the magnetization, thereby having a dominant effect on the magnetic response. The DW type studied here is the vortex-like DW in cylindrical nanotubes.

In the present paper, we investigate the dynamics of a DW in magnetic nanotubes under a magnetic field. We assume the magnetic manotubes to be confined within ideal cylindrical shapes, as shown in Fig. 1. First, we add a magnetic field along the z axis and establish the equation of motion of a domain wall by using the Landau–Lifshitz–Gilbert equation in magnetic nanotubes. In terms of a trial function we obtain the analytical solutions for domain-wall dynamics. Secondly, we get the domain-wall velocity, displacement, and distortion, which are shown as functions of time, respectively. At last, we discuss these quantities and energy variation, respectively. We find that the thickness of the nanotubes plays an important role in these quantities, such as the smaller thickness and the greater velocity of the domain wall.

Fig. 1. Vortex-like DW and cylindrical coordinates used in the calculation.
2. Dynamics of magnetization in magnetic nanotubes

The dynamics of magnetization is described by the Landau–Lifshitz–Gilbert equation[33,34] where γ is the gyromagnetic ratio, Ms is the saturation magnetization, and α is the Gilbert damping factor. Heff = −(δE/δm)/μ0 is the effective field with E being the system energy and μ0 being the vacuum permeability. For the magnetic nanotubes in Fig. 1 the local magnetization m takes the form , where θ represents the angle between the magnetization vector and the z axis, and φ is the out-of-plane angle of the magnetization vector projected in the φφ plane. The energy function[35] of magnetic nanotubes in Fig. 1 can be written as where A is the stiffness constant, K is the uniaxial anisotropy, β is the ratio of inner radius and outer radius, H is the applied external magnetic field along the z axis, and with l being the exchange length of material.

With the characteristic time t0 = μ0MsR2/(γA) and length , we can simplify Eq. (1) into the dimensionless form where heff = −(δε/δm) with the dimensionless energy ε = ∫ d z{(∂θ/∂z)2κ cos2 θ + η (1 + 2 cos2 φ) sin2 θh cos θ}, where κ = KR2/A, η = (2log1/β)/(1 − β2), and h = AH/(μ0R2Ms). Substituting the dimensionless energy ε into Eq. (3) we obtain where α1 = 1 + α2. With the help of Eqs. (4) and (5) we will investigate the magnetization motion in magnetic nanotubes driven by the external field.

3. Domain wall dynamics

The Walker solution analysis[3638] has been extensively adopted to investigate the moving DW in response to a magnetic field[3942] or spin-polarized current.[4350] The main idea of the so-called Walker method is to simplify the DW motion into a dynamic problem described by the center position and rotation angle. Therefore, the dynamics of a single DW can be qualitatively comprehended by one-dimensional analytical models that predict a rigid-body propagation below the Walker breakdown and an oscillatory DW motion above it.[3941] To this purpose, we follow Walker’s analysis of domain-wall motion by introducing a trial function where we have assumed that the projection of the magnetization vector in the domain wall on the ρφ plane only depends on the time. The second equation in Eq. (6) shows that the domain wall motion is characterized by the variation of width W(t) and the wall center with z0(t). That is to say, d z0(t)/dt stands for the domain wall velocity V(t), and z0(t) stands for the domain wall displacement.

Substituting this trial function in Eq. (6) into Eqs. (4) and (5), we get The above equations are the ordinary first-order differential equations, which determine the domain-wall width and the rotation of the domain-wall plane. Equations (7)–(9) can be solved by the four-order Runge–Kutta method and the numerical results are represented in Fig. 2, where the domain-wall velocity, displacement, and distortion are shown as functions of time.

Fig. 2. (color online) The DW angle ϕ, width W, velocity V, and the displacement z0 as a function of time t. The parameters are κ = 0.3, α = 0.1, and h = 0.5.

From Fig. 2(a), we observe that ϕ changes from the initial value gradually to the maximum, which is affected by the ratio of inner radius and outer radius β. When β is larger, i.e., the smaller the thickness of the nanotubes, the larger the maximum value of ϕ is. From Fig. 2(a) we find that dϕ/dt = 0 as t → ∞, and the domain wall distortion is maximum. Therefore, from Eq. (8) we obtain which implies the maximum reduction of the domain-wall width as It denotes that the domain-wall width decreases gradually, and the smaller the thickness of the nanotubes, the greater the width of the domain wall, i.e., the smaller the distortion of the domain wall structure when the domain wall structure is stable, as shown in Fig. 2(b). We can determine the nanotubes thickness to increase the stability of domain wall motion. The above results in Eq. (11) also show that in magnetic nanotubes the domain wall width is determined by the stiffness constant, anisotropy constant, and cylindrical symmetry (inner radius and outer radius).

In the initial time t = 0, ϕ = 0, and the initial DW velocity is zero. Under the conjunctive action of external field and damping the DW velocity attains the maximum at t → ∞ which is in linear proportion to the external field h. When β takes the different values, the DW velocity variation during the application of the external field is shown in Fig. 2(c). We find that when β is small, the DW velocity grows very quickly at first and the discrepancy for the different value β is very small. While the velocity with a larger β will gradually surpass that of the smaller β. When the domain wall structure is stable, the smaller the thickness of the nanotubes, the greater the velocity of domain wall. In other words, a thinner thickness of the nanotubes makes the spread of domain wall on nanotubes stronger. We also see that the DW structures in nanotubes are similar to those of flat strips, and have the significantly superior stability by reducing the thickness of the nanotubes. From Eq. (12) we can get the expression of the velocity when t → ∞, and the DW velocity is proportional to the width.

We estimate this displacement by integrating the velocity From Fig. 2(d), we observe that the DW displacement grows continuously with time linearity. Notice that the DW displacement is inversely proportional to the damping constant and the resistance of DW displacement implies the material permeability and coercive force. Domain wall displacement resistance mainly comes from the internal stress of material internal, doping, porosity, and defects on the domain wall pinning. Under the action of the magnetic field, the domain wall can move steadily. That is to say, the magnetic field can always drive domain wall motion.

As shown in Fig. 3(a), the domain wall width becomes broadened with the increasing ratio of inner radius and outer radius β. The domain wall width is determined by the balance between the exchange interaction and anisotropy, in which the exchange interaction favors the wide width of DW, while the anisotropy tends to make the DW width narrow. That is to say, the magnetization reversal needs a smaller external energy supplement in the magnetic nanotube. From Fig. 2(b), we can find the DW width increases with the increasing ratio of inner radius and outer radius β, while the velocity increases with the increasing ratio of inner radius and outer radius β, as shown in Fig. 3(b). This means we can accelerate domain wall motion like one-dimentional nanostructure and two-dimentional nanostructures, not only by increasing the magnetic field, but also by reducing the thickness of the nanotubes.

Fig. 3. (color online) The curve of domain wall width W and velocity V with the change of β. The parameters are κ = 0.3, α = 0.1, and h = 0.2.

It is interesting to study the energy variation under the field, i.e., dE/dt = Pα + Ph, with Pα = −h2W(∞)/α2 and Ph = 0. The first term Pα is the intrinsic damping power due to all kinds of damping mechanisms described by the phenomenological parameter. Ph is the external power. Note that the intrinsic damping power is always negative, whereas the total energy change rate is proportional to the negative DW velocity. Thus, energy is absorbed by the external field source. If we apply a constant magnetic field, the system is under the influence of the damping effect, the energy will gradually decrease. This shows that the damping effect is dissipated energy.

4. Conclusions

In this paper, we have solved the Landau–Lifshitz–Gilbert equation with the magnetic field. We propose a trial function to analytically obtain solutions for domain-wall dynamics. Finally, the distortion is maximum and the exact rigid-body solutions are obtained analytically. The DW can achieve greater velocity in magnetic nanotubes than the rigid-body DW in ferromagnetic nanowires driven by the same magnetic field. We find the velocity becomes bigger with the increasing ratio of inner radius and outer radius. That is to say, we can accelerate domain wall motion, not only by increasing the magnetic field, but also by reducing the thickness of the nanotubes.

Reference
[1] Parkin S S P Hayashi M Thomas L 2008 Science 320 209
[2] Allwood D A Xiong G Cooke M D Faulkner C C Atkinson D Vernier N Cowburn R P 2002 Science 296 2003
[3] Hayashi M Thomas L Moriya R Rettner C Parkin S S P 2008 Science 320 209
[4] Xiong G Allwood D A Faulkner C C Atkinson D Petit D Cowburn R P 2005 Science 309 1688
[5] Parkin S S P 2009 Scientific American 300 76
[6] Khizroev S Kryder M H Litvinov D Thompson D A 2002 Appl. Phys. Lett. 81 2256
[7] Hulteen J C Martin C R 1997 J. Mater. Chem. 7 1075
[8] Li F S Zhou D Wang T Wang Y Song L J Xu C T 2007 J. Appl. Phys. 101 6030
[9] Nielsch K Castano F Matthias S Lee W Ross C 2010 Adv. Eng. Mater. 7 217
[10] Sang J S Xia B Nan A Ghandehair H Sang B L 2006 J. Controlled Release 114 143
[11] Zhu Y Yang Q Zheng H G Gao L Yang Z Qian Y 2006 Mater. Chem. Phys. 96 506
[12] Berger L 1986 Phys. Rev. 33 1572
[13] Lee J Y Lee K S Choi S Guslienko K Y Kim S K 2007 Phys. Rev. 76 184408
[14] Slonczewski J S 1996 J. Magn. Magn. Mater. 159 L1
[15] Guslienko K Y Lee J Y Kim S K 2008 IEEE Trans. Magn. 44 3079
[16] Sun Z Z Schliemann J 2010 Phys. Rev. Lett. 104 037206
[17] Yan P Wang X R 2009 Phys. Rev. 80 214426
[18] Tretiakov O A Clarke D Chern G W Bazaliy Y B Tchernyshyov O 2008 Phys. Rev. Lett. 100 127204
[19] Beach GS Nistor C Knutson C Tsoi M Erskine J L 2005 Nat. Mater. 4 741
[20] Yang J Nistor C Beach G S D Erskine J L 2008 Phys. Rev. 77 014413
[21] Wang X R Yan P Lu J He C 2009 Ann. Phys. 324 1815
[22] Atkinson D Allwood D A Xiong G Cooke M D Faulkner C C Cowburn R P 2003 Nat. Mater. 2 85
[23] Wang X R Yan P Lu J 2009 Euro. Phys. Lett. 86 67001
[24] Kajiwara Y Harii K Takahashi S Ohe J Uchida K Mizuguchi M Umezawa H Kawai H Ando K Takanashi K Maekawa S Saitoh E 2010 Nature 464 262
[25] Chumak A V Serga A A Jungfleisch M B Neb R Bozhko D A Tiberkevich V S Hillebrands B 2012 Appl. Phys. Lett. 100 224403
[26] Tatara G Kohno H 2006 Phys. Rev. Lett. 96 086601
[27] Zhang S Li Z 2004 Phys. Rev. Lett. 93 127204
[28] Schryer N L Walker L R 1974 J. Appl. Phys. 45 5406
[29] Otaora J A Lopez-Lopez J A Núñez A S Landeros P 2012 J. Phys.: Condens. Matter 24 436007
[30] Otaora J A Lopez-Lopez J A Vargas P Landeros P 2012 Appl. Phys. Lett. 24 436007
[31] Yan M Andreas C Kakay A Garcia-Sanchez F Hertel R 2011 Appl. Phys. Lett. 99 190
[32] Kakay A Yan M Andreas C Hertel R 2013 Phys. Rev. 88 220412
[33] Gilbert T L 2004 IEEE Tran. Magn. 40 3443
[34] Skomski R Coey J M D 1999 High Magnetic Fields New York
[35] Landeros P Á S 2010 J. Appl. Phys. 108 033917
[36] Schryer N L Walker L R 1974 J. Appl. Phys. 45 5406
[37] Slonczewski J C Malozemoff A P 1979 Academic New York 309
[38] Thiaville A Nakatani Y 2006 Berlin Springer III 161
[39] Beach G S Nistor C Knutson C Tsoi M Erskine J L 2005 Nat. Mater. 4 741
[40] Lu J Wang X R Yan P He C 2009 Ann. Phys. 324 1815
[41] Wang X Yan P Lu J 2009 Europhys. Lett. 86 67001
[42] Sun Z Z Schliemann J 2010 Phys. Rev. Lett. 104 037206
[43] Tsoi M Fontana R E Parkin S S P 2003 Appl. Phys. Lett. 83 2617
[44] Tatara G Kohno H 2006 Phys. Rev. Lett. 96 086601
[45] Zhang S Li Z 2004 Phys. Rev. Lett. 92 207203
[46] Thiaville A Nakatani Y Miltat J Suzuki Y 2005 Europhys. Lett. 69 990
[47] Saitoh E Miyajima H Yamaoka T Tatara G 2004 Nature 432 203
[48] Tatara G Saitoh E Ichimura M Kohno H 2005 Appl. Phys. Lett. 86 4825
[49] Zhang S Li Z 2004 Phys. Rev. Lett. 93 127204
[50] He P B Xie X C Liu W M 2005 Phys. Rev. 72 172411